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ABSTRACT 

The motion of a point like object of mass M passing through the background potential of 
massive collisionless particles (m << M) suffers a steady deceleration named dynamic friction. 
In his classical work, Chandrasekhar assumed a Maxwellian velocity distribution in the halo and 
neglected the self gravity of the wake induced by the gravitational focusing of the mass M. In this 
paper, by relaxing the validity of the Maxwellian distribution due to the presence of long range 
forces, we derive an analytical formula for the dynamic friction in the context of the q-nonextensive 
kinetic theory. In the extensive limiting case {q = 1), the classical Gaussian Chandrasekhar result 
is recovered. As an application, the dynamic friction timescalc for Globular Clusters spiraling to 
the galactic center is explicitly obtained. Our results suggest that the problem concerning the 
large timescale as derived by numerical A^-body simulations or semi- analytical models can be 
understood as a departure from the standard extensive Maxwellian regime as measured by the 
Tsallis nonextensive g-parameter. 

Subject headings: Dynamical Friction, Nonextensivity, Globular Clusters 



1. Introduction 

It is widely known that a massive object of mass 
M such as a Globular Cluster passing through 
a background of non-colliding particles suffers a 
gravitational effect usually referred to as Dynam- 
ical Friction (DF). Historically, this problem was 
first studied by Chandrasekhar (1943) who first 
recognized its true dissipative nature. He ana- 
lyzed the idealized case where a point mass moves 
through an infinity, homogeneous sea of field par- 
ticles and showed that a fraction of the kinetic 
energy of the incoming object is transferred to the 
stellar collisionless population whose distribution 
was described by a Maxwelllian velocity. 

The DF mechanism is now a classical ef- 
fect for description and evolution of almost all 
many-body astrophysical systems. Some exam- 
ples involve the formation of stellar galactic nu- 
clei via merging of old Globular Clusters (GCs) 
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([Tremaine et al. 19751) , the transformation from 
non-nucleated d warf galaxies into nucleated ones 
(|Oh fc Lin 200n[ ). :he behavio r of radio galaxies 
in galaxy clusters (Nath 20081 ) . nonlinear gaseous 
medium ( Kin Sz Kin 2009j) and the fi eld particles 
with a mass spectrum ( Ciotti 20101 ). Tradition- 
ally, such investigations were carried out in the 
framework of Newtonian gravity, however, alter- 
native gravity theories like the Modified Newto- 
nian Dynamics (MO ND) has also been considered 
(|Nipoti et al. 20081 ). 

In the last few years, several authors have dis- 
cussed the problem related to the dynamical fric- 
tion timescale (tdf) of a GC orbiting dwarf galaxies 
or infalling satellite galaxies in clusters (Read et al. 
2006; Goerdt et al. 2006; Sanchez-Salsedo et al. 
2006; Nath 2008; Cowsik et al. 2009; Inoue 2009; 
Namouni 2010, Can et al. 2010). In particular, 
when dwarf galaxies have a cored dark matter halo 
with constant density distribution in its center it 
has been found that the DF effects will be consid- 
erably modi fied. Some ana lyses b ased on TV- body 
simulations ( Goerdt et al. 2006; I noue 20091 ) and 
semi-analytical models have shown that the ex- 
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pected sinking of CGs to the galactic center may 
take a time beyond the age of the universe. 

In order to explain this suppression effect 
and the consequent problem of extremely long 
tdf, some authors introduced additional ingre- 
dients, among them: the loss mass of the GCs 
by stripping and heating t idal ( Gan et al. 20ld : 
Bovlan-Kolchin et al. 20081 ) . and the interac- 
tion between th e dark matter halo and the GC 
(|0h fc Lin 2OOOI: iRead et al. 20061: llnoue 20091 ). 
All these possibilities are based on the standard 
Chandrasekhar's DF formula derived by assuming 
a Maxwellian distribution. Here we investigate a 
different route by considering that the solution for 
this problem is related to a proper extension of 
the underlying statistical approach. 

It is well known that the so-called nonextensive 
statistical approach provides an analytical exten- 
sion of Boltzmann-Gibbs (BG) statistical mechan- 
ics which is very convenient when long-range forces 
are present and/or the system is out (but close) to 
a thermal equilibrium state. This ensemble the- 
ory is based on the formulation of a generalized 
entropy proposed by TsalHs (1988,2009) 



Sq = ks 



1 V^VK i 



(1) 

which reduces in the limit g — !■ 1 to the BG entropy 
Sbg — X^ll^i-P* In Pi, since pi is the probabil- 
ity of finding the systems in the microstate i, W 
is the number of microstates and fcs is the Boltz- 
mann constant. However, when the index 5 7^ 1, 
the entropy of the system is nonextensive, i.e, 
given two subsystems A and B, the entropy is 
no more additive in the sense that Sq{A + B) = 

SqiA) + SqiB) + (1 - q)SqiA)SqiB). Thc loUg- 

range interactions are associated to the last term 
on the r.h.s. which accounts for correlations be- 
tween the subsystems with the index q quantifying 
the degree of statistical correlations. Such a sta- 
tistical description has been successfully applied 
to many complex physical systems ranging from 
physics to astrophysics and plasma physics, among 
them: the electrostatic plane- wave propagation in 
a collisionless thermal plasma (Lima, Silva & San- 
tos 2000), th e peculiar velocity fu nction of galax- 
ies clusters ( Lavagno et al. 19981). g ravothermal 
instability ( Taruva Sz Sak agami 20021 ) . the kinetic 
concept of Jeans gravitational instability (Lima, 
Silva & Santos 2002), and the radial and projected 



density profiles for t wo large classes of iso ther- 
mal stellar systems (jLima fc de Souza 20051 ). A 
wide range of physical applications can also 
be seen in Gell-Mann & Tsallis 2004 (see also 
'http : / / tsallis . cat . cbpf . br /biblio.htm| for an up- 
dated bibliography). 

In this letter, by assuming that a self-gravitating 
collisionless gas is described by the nonextensive 
kinetic theory (Silva et al. 1998; Lima et al. 2001), 
we derive a new analytical formula for dynamical 
friction which generalizes the Chandrasekhar re- 
sult. As an application, the DF timescale (tdf) for 
GCs falling in the galaxies center is derived for 
the case of a singular isothermal sphere. This re- 
sult suggest that the long timescales for GCs can 
be understood as a departure from the extensive 
regime. In other words, there is no suppression of 
the DF since the long time can be just the sta- 
tistical price to pay by the presence of long range 
forces acting on the gravitational systems. 

2. Dynamical Friction and Nonextensive 
Effects 

By following Chandrasekhar (1943), the DF de- 
celeration on a test mass M moving with velocity 
vm in a homogeneous and isotropic distribution 
of identical field particles of mass m and number 
density no reads: 



M 



dt 
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VM, 



(2) 

where G is the gravitational constant, m is 
the mean mass of field stars and f{v) repre- 
sents their velocity distribution. The parameter 
A — Pmax/Pmin depends on the ratio of the the 
maximum (pmax) and minimum (pmin) impact 
parameters of the encounters contributing to gen- 
erate the dragging force. 

In the applications of DF, it is usually assumed 
that the distribution function of the stellar veloc- 
ity fiel d can be described by a Maxwellian distri - 
bution (jBinnev fc Tremaine 20"08HFellhauer 20081 ) 



fix.) = 



HQ 
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(3) 



where AT^ = v/\/2a denotes a normalized velocity 
with a indicating their dispersion. The integration 
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of (HI) results: 

dvM _ AttIu KG'^Mp{r) 



dt 



Hi{Xm)vm, (4) 



where p{r) — n^m and the function Hi(Xm) is 
given by 



Hl(XM) = er/(XM)-^e-^^ 



with erf{XM) defining the error function as 
2 rXM 



erfiX 



e-^'*dX^. 



(5) 



(6) 



Now, in order to investigate the nonextensive 
effects on the Chandrasekhar theory, let us con- 
sider that the stellar field obeys the following 
power-law (Silva, Plastino & Lima 1998, Lima, 
Silva & Plastino 2001, Lima & de Souza 2005): 
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where the so-called (7-exponential is defined by 



(7) 



(8) 



and the quantity Aq denotes a normalization con- 
stant which depends on the interval of the q- 
parameter. For values of q < 1, the positivi- 
ness of power argument means that distribution 
above exhibits a cut-off in the maximal allowed 
velocities. In this case, all velocities lie on the 
interval {0,Vmax) and their maximum value is 
Vmax = v^cr/Vl ~ Q- Taking this into account 
one may show that the normalization constant Aq 
can be written in terms of Gamma functions as 
follows: 



^ 1 — q ' 

A, = (g-l)3/2pP#l4-, q>l 



(9) 



For generic values of (? ^ 1, the DF 
is a power law, whereas for g = 1 it reduces 
to the standard Maxwell-Boltzmann distribu- 
tion function ([3]) since — )• 1 at this limit. 
Formally, this result follows directly from the 
known identity, limd-»o(l -|- dy)d = exp(y) 
( Abramowitz fc Stegun 19721 ). The distribution 
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Fig. 1.— The Hq{X) function. The sohd black 
curve is the result based on Chandrasekhar the- 
ory (Hi{x)). The remaining curves show the q- 
corrections for several values of the q-index. 

([7]) is uniquely determined from two simple re- 
quirements (Silva et al. 1998): (i) isotropy of 
the velocity space, and (ii) a suitable nonexten- 
sive generalization of the Maxwell factorizability 
condition, or equivalently, the assumption that 
/(^) 7^ /(^a;)/(?^i/)/(f2:). The kiuetic foundations 
of the above distribution were also investigated 
in a deeper level through the generalized Boltz- 
mann's equation, in particular, it was also shown 
that the kinetic version of the Tsallis entropy satis- 
fies an extended TJ^-theorem (Lima, Silva & Plas- 
tino 2001). 

Now, by considering that the power-law dis- 
tribution ([Tj is a valid description for the stellar 
velocity distribution we conclude that the expres- 
sion describing the DF in this extended framework 
takes the following form: 



M 



167r2G'2(lnA)Mp(r) 



dt 



Aq X 



/ X^eq{X,)dX^VM, (10) 
"'0 



which after an elementary integration can be 
rewritten as: 

dvM 4^G2lnAp(r)M f 2 , „ ,^ , 



dt 



b-iq 



(11) 



where Hq{XM) is the general function depending 
on the q-parameter (compare with Eq. (5)) 



Hq{XM) = I,{Xm) - ^^Aqe^-^XM) 



(12) 



3 



In the above expression, the integral 



2A„ 



IqiXu) = -j^ / eq{X,)dX,, (13) 



is the g-generahzation of the error function (see 
Eq.(6)). 

As one may check, the nonextensive expression 
for the DF (including the auxiliary functions Hq 
and Iq) reduces to the Chandrasekhar result in the 
Gaussian limit (g — >■ 1). It shows clearly that the 
collective effect from gravitational interactions of 
M (with all stars of the field) is strongly depen- 
dent on the statistical model. An interesting as- 
pect of the above formulae is that the results are 
expressed by analytical expressions. In principle, 
they can be useful for semi- analytical implementa- 
tions because the easy comparison with the stan- 
dard approach (see next section). Naturally, we 
are also advocating here that the idealized frame- 
work based on the Maxwellian distribution (Chan- 
drasekhar 1943) may be in the root of some the- 
oretical difficulties shown by iV-body simulations, 
like the ones related to the decay orbits of GCs. 

3. Decay of Globular Orbits 

In order to illustrate some consequences of the 
above derivation, let us now analyze the nonexten- 
sive solution for the decaying orbit of a GC in the 
stellar galactic field. As a GC orbits through the 
galaxy field, it is subject to DF due to its inter- 
action with the stellar distribution. By assuming 
spherically symmetric star distribution, the drag- 
ging force decelerates the cluster motion which 
loses energy thereby spiraling toward the galaxy 
center. Therefore, whether the GC is initially on a 
circular orbit of radius , it is convenient to define 
an average DF timescale, t^/, as the time required 
for the cluster reach the galaxy center. For the 
sake of simplicity, we also consider that the mass 
density distribution of the galaxy is described by 
the singular isothermal sphere 



p(r) = 



1 



47rG' V r 



(14) 



with Vc being circular speed and a = Vc/V^ the 
velocity dispersion. This simplified mass distri- 
bution has the benefit of having a planar rotation 
curve and therefore might be considered as a crude 



but minimally realistic distribution for the exter- 
nal region of normal galaxies. Following standard 
lines, the frictional force undergone by the clus- 
ter of mass M moving with speed Vc through the 
stellar field now reads: 



F : 



2 \ /Af 

^GlnA( — ) Hq{l), (15) 
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where Hq{\) is the general function ([T^ written in 
the coordinate X — {vc/cFy/2) — 1. Note also that 
the integral Iq{XM) defined in ([T^ now reduces 
to 



, , 2A„ / 1 1 3 



(16) 



where 2Fi{a,b;c] z) is the Gauss hypergeometric 
function. Either from the above representation 
or from the integral form (|13|) . we see that the 
error function er/(l) is obtained as a particu- 
lar case in the e xtensive regime, that is, /i (1 ) = 
er/(l) w 0.8427 (jAbramowitz fc Stegun 19721 ). It 
means that H^{l) = er/(l ) - (2/v^)e-i « 0.428 
(jBinnev fc Tremaine 20081 ) . 

Now, returning to expression ()15p . we recall 
that the dragging force is tangential to the clus- 
ter orbits, and, therefore, the cluster gradually 
loses angular moment per unit mass L at a rate 
dL/dt ~ Fr/A'I. Since L = rvc we can rewritten 
equation (fT5|) as 



dr 
'dt 



5-3q 



— )lnAff,(l). (17) 

Vc J 



By solving this differential equation subjected to 
the initial condition, r(0) — ri, we find that the 
cluster reaches the galaxy center after a time 



t 



(9) _ 
df - 



5-3(7 



GMlnAHq{l)' 



(18) 



This nonextensive timescale for decaying orbits 
of GCs generalizes the Chandrasekhar result (see 
Binney fc Tremaine 2008) which is readily recov- 
ered in the Gaussian extensive limit {q — 1). 

At this point, it is interesting to compare the 
above nonextensive prediction with the standard 
result based on the Chandrasekhar approach. To 
begin with, let us assume typical values for the 
parameters r^, Vc and M, namely: Tj — 2Kpc, 



= 250fcr 



and M = lO^'Mtr,. With these 
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Fig. 2. — Behavior of the relative time scale ra- 
tio Tq. We see that for q > 1, the characteris- 
tic nonextcnsivc time scale for dynamic friction 
can be much greater than in the standard Chan- 
drasekhar approach. 

choices we get: 

(q) _ l.Mx lQii ( 5-3q \ f n y 
~ Hg{l)lnA \ 2 J \2kpcJ 

which reduces to the standard value in limiting 
case {q = 1) as given by Binney & Tremaine 
(2008). The nonextcnsivc corrections arc more di- 
rectly quantified by introducing the dynamic time 
ratio, T{q) = t^^f /t'^Jj} , where t^^^ denotes the 
Chandrasekhar result. By using (18) we find 

^<') K^) til 

where the function Hg{X) was defined by Eq. 
(11). 

In Figure 1, we display the nonextcnsivc cor- 
rections for a large range of the nonextcnsivc 
g-parameter. As a general result, we see that 
the T(q) ratio is strongly dependent on the q- 
parameter. The nonextensive time scale is greater 
or less than the extensive Chandrashckhar re- 
sult depending on the interval of the q-parameter. 
Note also that t'^ is greater or smaller than t'^Jj} 
if 9 > 1 or g' < 1, respectively. 

4. Conclusions 

We have derived the g-dynamic friction force 
for a point mass moving through a homogeneous 



background in the context of the nonextcnsivc ki- 
netic theory. Simple and analytical forms were ob- 
tained, and, as should be expected, they smoothly 
reduce to the standard Chandrashckar results in 
the extensive limiting case (g = 1). However, 
for g 7^ la large variety of qualitatively dif- 
ferent behaviors arc predicted when the free pa- 
rameter g is continuously varied (see Figs. 1 e 
2). As an application, we have discussed the dy- 
namical timescalc for a glob\.ilar cluster collapsing 
to the center of a massive dark matter halo de- 
scribed by an isothermal sphere. The results pre- 
sented here suggest that the problem related to the 
large timescalc shown by numerical iV-body simu- 
lations and semi-analytical models may naturally 
be solved (with no ad hoc mechanism) by taking 
a proper g-nonextensive distribution with param- 
eter greater than unity. Applications to more re- 
alistic density profiles like the lowered nonexten- 
sive halos distribution (Silva, de Souza & Lima 
2009; Cardone, Leubner & Del Popolo) and a de- 
tailed comparison with semi-analytical models will 
be discussed in a forthcoming communication. 
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